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Abstract. - We analyze the morphological transition of a one- dimensional system described 
by a scalar field, where a flat state looses its stability. This scalar field may for example account 
for the position of a crystal growth front, an order parameter, or a concentration profile. We 
show that two types of dynamics occur around the transition: weakly nonlinear dynamics, or 
highly nonlinear dynamics. The conditions under which highly nonlinear evolution equations 
appear are determined, and their generic form is derived. Finally, examples are discussed. 



In the study of pattern formation, weakly nonlinear equations play a central role. By 
construction, these equations catch the main effects of nonlinearities via a limited number of 
nonlinear terms added to a linear equation. This approach has been used for a wide variety of 
physical systems such as crystal growth [1], reaction-diffusion systems [2], flame fronts [3], or 
phase separation [4]. Weakly nonlinear equations can be derived from a multi-scale analysis 
when separation of scales is possible. This is for example the case in the vicinity of an 
instability threshold, where the system is weakly unstable, or in the analysis of amplitude 
and phase dynamics of modulated structures [5]. These equations are also obtained from 
renormalization techniques [6]. 

Some analysis and attempt of classification of generic nonlinear equations based on sym- 
metry or geometry have already been reported in the literature [5,7,8]. The most systematic 
approach up to now was that of Ref. [8] , where nonlinear equations result from the expansion in 
Cartesian coordinates of dynamics expressed in intrinsic coordinates. We here present a more 
general approach based on a multi-scale analysis. For the sake of simplicity, we assume that 
dynamics is local and that an instability appears at long wavelength at the instability thresh- 
old. From the assumption that the stabilizing or nonlinear terms do not scale with the small 
parameter of the expansion e, we find that the Benney, and the sand ripple [9] equations are 
expected in systems with translational invariance, and that the convective Ginzburg-Landau 
equation is expected in absence of translational invariance. 

Furthermore, our approach determines the range of validity of weakly nonlinear expansions 
even when e is present in the stabilizing or nonlinear terms. As a central result, we show that 
the weakly nonlinear approach breaks down for a large class of front dynamics. The main 
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property of these systems is to be in the vicinity of a variational steady-state (which may be 
thermodynamic equiUbrium). We obtain highly non- linear evolution equations in two cases: 
(1) Conserved dynamics with translational invariance. We provide explicit illustrations of 
this case in the context of Molecular Beam Epitaxy. (2) Dynamics without translational 
invariance, but with a translationally invariant stabilization process (such as surface tension). 
This situation for example applies to phase separation [4], or to amplitude equations for 
modulated patterns [5]. 

Unexpectedly, highly non- linear equations exhibit the "central" symmetry in case (1), and 
the left-right symmetry in case (2), which are not present in the original physical problem. 
Furthermore, a Lyapunov functional may be found in case (2), although we are in a fully 
non-equilibrium situation. 

Let us first consider the case where the front obeys translational invariance. This means 
that df h does not change when the whole front is translated via h — > h + Hq, where is a 
constant. Then, dynamics does not depend on the position of the front h{x,t), but only on 
its derivatives with respect to time t and space x. For a linear analysis, we consider a small 
perturbation corresponding to a unique Fourier mode h(x, t) = h^^jk exp(wi -I- ikx). Inserting 
this relation in any front dynamics model then provides w as a specific function of k. Assuming 
that dynamics is local, Ke[a;] and 3TO[ci;] respectively only involve even and odd powers of k. 
A long wavelength (i.e. small fc) expansion then leads to: 

^e[uj] = L2k^ + Lik'^ + o{k^); '^m[uj] = L^k^ + o{k^) (1) 

From translational invariance, the mode fc = is marginally stable, and thus, there is no 
constant term in Ke[a;]. Moreover, we have performed a Galilean transform x ^ x + Lit in 
order to eliminate the linear term Lik in 3m[a;]. As fc — ^ 0, one has L2k^ ^ L4k*. Therefore, 
the criterion for an instability to occur (i.e. ^e[uj] > 0) at long wavelength is simply L2 > 0. 
Since we restrict the analysis to instabilities occurring at long wavelength at the threshold, 
one should generically require L4 < 0. 

We now perform a multi-scale analysis in the vicinity of the instability threshold based on 
an expansion with the small parameter e L2 (^). From Ec^.l^, the unstable modes are those 
for which < k < kc = (— L2/L4)^/^ e^^^, and the most unstable mode is km = fcc/2^/^. 
Therefore, the unstable modes have k ~ e^/^ and ^e[uj] ^ e^, so that the relevant spatio- 
temporal scales are x ~ k~^ ~ and t ~ 5fte[a;]~^ ~ (^). Using Eq.QJ for the linear 
part, the general form of a weakly nonlinear equation is defined as: 

dth = -L2d,^h - Lsd^.^h + LA.x.h + t^[d^ndt]'[hr (2) 

In the last term, the brackets mean that we account at the same time for all terms containing 
n spatial and I temporal derivatives, and m times h (where m > 1), with arbitrary numerical 
prefactors. Translational invariance imposes that n + I > m. Moreover, we do not know the 
value of 7 a priori, which is a consequence of the specific properties of the system. 

We first analyse the self consistency of Eq.Q, using power counting arguments. Indeed, 
nonlinear terms and linear terms relevant to stability should be of the same order in e. Stating 
that /i ~ e", we then find: 

a= [2-7-n/2-2/]/(m-l). (3) 

(^)We could assume that L4 ~ e*", as long as the instability occurs at long wavelength (i.e. 64 < 1). Defining 
t' = e^it, 7' = [7 + &4{l — 1)]/(1 — 154), and e' = we then have again an equation of the form Eg. 01. 

(^) Although we do not present it here for the sake of clarity, we should introduce a propagative time-scale 
~ €~^^^ L^^ , related to the dispersive term. This would not change the results. 
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The weakly nonlinear equation |2Jl is the result of an expansion of dth -which is an unknown 
function of the derivatives of h~ in the limit e — > 0. In order to perform this expansion, 
we need all the derivatives of h to be smaller than one. Since ~ e" and x ~ e~^^^, we 
need pnt+n:^/2+a ^ i \ sufficient condition is that the largest derivative dxh is 

smaller than one. This reads: d^h - e"+i/2 ^ i, We shall thus require that a > —1/2. This 
condition allows the usual gradient expansion. Combining the different constraints, we find 
the condition: 

27 + 3(;- I) < m-n-Z < 0, (4) 

which determines all possible terms which may intervene in a weakly nonlinear equation. 

We have now determined the conditions I0J under which a weakly nonlinear expansion is 
well defined. We shall now deal with the more subtle question of the self-consistency of the 
dynamics resulting from a weakly nonlinear equation. The central question is the interplay 
between dominant and subdominant terms. We define the rescaled variables h = e~°'h, 
X = e^^^x, and t = e^t. Let us then consider two nonlinear terms Ni and N2 (with their 
rescaled forms iVi and N2) and the related values of a from Eq.||21l such that ai > 0,2 (■^)- 
Assuming that a = ai, one finds that 

Ar2/iVi ^ e(™^-i)("i-"2)(7V2/iVi) (5) 

-we recall that > 1. Depending on the dynamics in presence of iVi in Eq.Q, three cases 
may be observed: (i) The solution h(x, t) (or its derivatives for a system with translational 
invariance) is bounded. Then Ni are bounded. It follows that N2 ^ Ni. (ii) h(x,t) scales 
with time, i.e. the typical scales of the patterns evolve ash^ P'* and x ^ t°"^ . Then Ni ~ t'*', 
and N2/N1 - £(™2-i)(ai-a2)iK2-«;i^ Therefore, if ki > K2, we have iVi > N2 at all times. 
If Ki < K2, then N2 ^ Ni after a time tc ^^£-("i2-i)(ai-a2)^ where is the time for 
appearance of the instability from the linear analysis. As e — > 0, one finds that tc ^ t^. (iii) 
h{x,t) exhibits an exponential -or even faster- increase of the amplitude (including the case 
of singularities in finite time). Then N2 may become of the same order as A^i after a time 
tc ^ tL (with possible logarithmic corrections). 

From an inspection of the three cases mentioned above, we see that the only self-consistent 
choice in order to describe the dynamics at timescales that are much longer than is to find, 
if it exists, the nonlinearity with the biggest value of a. The dynamics may then be: (i) 
bounded, and the time tc under which the weakly nonlinear equation is valid is infinite; (ii) 
power-law in time. The time tc is then either infinite, or increasing as e to some negative 
power when e 0. But in the case (iii), where the growth of the amplitude is exponential 
or faster, the weakly nonlinear equation does not provide a satisfactory description of the 
dynamics in the nonlinear regime. 

We conclude that a nonlinear analysis in the regime where e is small requires an expansion 
of the form |(5J, but also relies on the knowledge of the dynamics of ^ with the dominant 
nonlinearities. There is to our knowledge no general analytical tool which could systematically 
determine whether the nonlinear dynamics belongs to one of these 3 classes of dynamics. 
Therefore, a numerical solution is in general needed. 

Let us now consider some precise examples. In general, the dominant contribution depends 
on the precise values of a, which themselves depends on 7. Therefore, the specific physical 
ingredients of the system will determine the most relevant terms. Nevertheless, we shall 
first consider the simplified case where 7 = 0. Relations then readily show that I = 0. 

{^)Ni, with 1 = 1,2 more generally describe a sum of nonlinear terms having the same value of ai. 
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From Eq.Q, we are now able list all possible nonlinearities in a weakly nonlinear expansion: 
{dxh)"^]dx{dxh)"^]dxx{dxhY^] {dxh)"^~'^{dxxh)^ , where m > 1. Using this result we find that 
the expected equation is the Benney equation: 

dth = -dxxh + Pdxxxh - dxxxxh + [dxhf (6) 

because it leads to the largest possible value of a, which is a = 1, and because the solution 
of Eq.® leads to a saturation of the amplitude. The variables x, t, h have been normalized 
in Eq.® so that only one constant remains. Eq.® is non- variational and exhibits order 
or chaos when the parameter (3 is larger or smaller than one respectively. In the chaotic 
limit (3 — 0, it is called the Kuramoto-Sivashinsky equation. The Benney and Kuramoto- 
Sivashinsky equations have been derived from multi-scale analysis in many physical situations, 
such as flame fronts [3], crystal step meandering [1] and bunching [10], or ion sputtering [11]. 

In Molecular Beam Epitaxy (MBE), atoms land on the surface but can usually not go back 
to the atmosphere. One then has from mass conservation: 

dth = F- dxj (7) 

where F is the incoming flux, and j is a mass flux along the surface. The constant term is 
eliminated by the transformation h h + Ft. Many other systems obey a conservation law, 
such as wave dynamics in thin liquid layers [12], or sand ripple formation [9]. 

Once again, we shall first assume that 7 = 0. From Eqs.O and the only nonlinear 
terms which are allowed are: dx{dxh)™'; dxx{dxh)™ . Following the same line as above, we 
find that the first dominant nonlinear term is dx{dxh)^. But this term leads to an increase 
of the amplitude which is faster than power-law [9]. This nonlinearity can be absent if it is 
proportional to , with 7 > 1/2, but also if some symmetry (such as h —h, or x ~x) is 
imposed. We then find: 

(8) 

which corresponds to a = 0. In Eq.®, t, x, and h have been normalized, and q are constants. 
This equation was first obtained by Csahok et al [9]. When f3 = and C2 = (e.g. when 
dynamics is variational, or when it exhibits the {h,x) — > {—h,—x) symmetry), and when 
C3 < 0, one recovers the Cahn-Hilliard [4] equation with a double well potential, which is 
known to lead to logarithmic coarsening [13]. In this case, the amplitude dxh remains finite 
at all times. Therefore, the expansion is self consistent, and higher order terms are negligible. 
As shown in Ref. [9], this behavior seems to persist when (3 and C2 are both non-zero. But 
when €3 = and C2 ^ 0, power law coarsening is found, with the wavelength t^l"^ and 
the amplitude ~ t^l"^ [14]. When C3 > 0, the dynamics leads to a local blow up of the slope 
dxh. Eq.||SI) describes a wide variety of systems, such as phase separation [4], sand ripple 
formation [9], and step bunching [14]. 

Up to this point, we have found that weakly nonlinear expansions are generically obtained 
from multi-scale analysis. This result was based on the assumption that 7 = 0. Nevertheless, 
when 7^0, the set of nonlinear terms which must be kept can change. For example, in 
Ref. [14], the step bunching instability on a vicinal surface is studied under growth and 
mobile atom migration. In this case, e is proportional to the growth rate for some given 
value of the migration rate. The terms dx{dxh)^ and dx{dxh)^ are absent in Eq.©, because 
they are multiplied by e'^ with 7 = 1. A more drastic consequence of a non- vanishing 7 is 
the possible break-down of the weakly nonlinear expansion. In the following, we indeed show 



dth = dx 



'dxh + Pdxxh - dx {dxhf + C2dx{dxhf 
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that conserved dynamics, in the vicinity of thermodynamic equihbrium, forbids any weakly 
nonhnear expansion. 

Let us first consider the conserved modeh dth = dx[Mdx{ST/6h)] , where 6 denotes the 
functional derivative, and J- = J dxcj). M > and (f> are function of the spatial derivatives of 
h. J-', which plays the role of an energy, and is a Lyapunov functional (i.e dtJ-^ < 0). Assuming 
that J-' is minimum for a straight front {dxh = 0), any initial perturbation would decay. This 
model may for example account for relaxation towards thermodynamic equilibrium. When it 
is driven by a small non-equilibrium force / > (not breaking mass conservation), the system 
will respond by an additional flux J. To leading order in /, the new dynamics reads: 

+oin, (9) 

where J is a function of the spatial and temporal derivatives of h. Linearizing Eq.Q, we find: 

dth = fJidxxh + fJ2dxxxh - (Mq - fJz)dxxxxh + ... (10) 

where Ji = do^hJ, J2 = do^^hJ, J3 = do^^^uJ , and Mq = Mdo^ha^h4' in the steady-state 
configuration. From the analogy between Eas. (|l()|l and QJ, we conclude that an instability 
occurs if Ji < 0, and that e ~ /. From the last term of Eg. 1)10(1 one then finds that L4 
does not scale with e. We now look for possible weakly nonlinear terms. For the first term 
in the brackets of Eq.©, Z = 0, 7 = 0, and n > m -|- 3, which is in contradiction with Q). 
If a term comes from the non-equilibrium contribution J, one has 7 = 1, and n > m + 1, 
which is again in contradiction with Finally higher order terms have 7 > 2, which is also 
in contradiction with Therefore, no nonlinear term satisfies (0J, and the small gradient 
constraint {dxh <^ 1) has to be waived. Thus, we choose a = —1/2, and since dxh ~ 1, the 
full nonlinear dependence of M, ST/6h, and J on dxh must be kept, while terms such as dxxh 
or dth are negligible. Therefore, J ^ A, M B, and ST/Sh — > dxC, where A, B and C are 
functions of dxh only. To leading order in e, Eq.® then takes the highly nonlinear form: (^) 

dth = dx[BdxxC + eA] (11) 

Unexpectedly, Eg. 1(11(1 is invariant under the "central" symmetry {h^x) {—h, —x), although 
the starting point Eq.Q is not. We shall notice that a formal expansion of A, B, and C with 
respect to dxh in Eg. 1(11(1 leads to an equation of the form (O with an infinite number of 
nonlinear terms, which all have the same value of a. Therefore, highly nonlinear equations 
can be considered as a special case of and the classification of the dynamics into the 3 
classes, bounded, power-law, and exponential, is still valid. 

An equation of the form 1(11(1 was first derived from a multi-scale analysis of crystal step 
meandering during MBE [15]. Although Ea. ((ll(l is not variational (except in some special 
cases, e.g. when Bdxh — A oy C — B) ^ the simple structure of its steady-states allows one to 
analyze in details its coarsening dynamics [19] as found in the case of step meandering from 
several recent works [15-17]. Eg. l(ll() was also introduced as a model for mound formation 
during MBE [18]. Our analysis now provides a frame to understand its origin. 

We now turn to the dynamics of a front without translational invariance. This situation 
occurs when a front, such as the free surface of an thin adsorbate, is subject to an external field 
which is a function of h (due for example to the substrate). It may also account for unstable 
concentration profiles in reaction diffusion systems [2], or for phase separating systems [4]. 

('*)A more general form for conserved dynamics when a = —1/2 is dth = dx[B3dxxm + 32(8x1^1)^ + 318x^1 + 
eBom], where Bi are function of m = dxh such that Bi — > constant when m — > 0. 



dth = dx 



Mdx^-^ + fJ 
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Since we use a similar approach to that presented above, we will be more concise here. We 
shall assume that there exist a flat steady state h — 0. At long wavelength, 

sReM « Lq + L2P. (12) 

The instability now occurs if Lq > 0, and we expect L2 < 0. We therefore choose e ~ Lq and 
the relevant spatio-temporal scales are: x e~^/^ and t ^ e^^. As in the previous case, the 
term Lik, which appears in Smfo;], can be cancelled with the help of a Galilean transform. 
The contribution 9m [w] ~ Lafc^ ~ e^/^ is negligible when L3 does not scale with a negative 
power of e. The general form of a weakly nonlinear equation is now: 

dth = Loh - L^d^^h + e^[d,r[dtf[hr (13) 

with TO > 1. The nonlinear term will be of the same order as the linear terms if: a = 
(1 — n/2 — 7 — O/l'^ ~ !)■ The condition for having weakly nonlinear dynamics is /i <C 1, 
which implies that a > 0. This leads to some restriction, namely: 

7<l-n/2-L (14) 

Once again, we assume that 7 = 0. Then I — 0, and the only possible nonlinear terms are: 
h™; dx{h™). The generic first contribution is h? . It leads to the Fisher-Kolmogorov equation, 
which has been extensively used for its traveling wave solutions [20]. But here, we must start 
from a front with zero average height, in which case h locally diverges in finite time. As before, 
the presence of a pref actor e"*", with 7 > 1/2, or the existence of symmetries such as h —> —h 
or {h,x) ■—>■ {~h, —x) may forbid this term. One then finds the convective Ginzburg-Landau 
equation 

dth = h + dxxh + (jh^ + fihdxh (15) 

where x h, and t are normalized. When a > 0, h again locally diverges in finite time. We 
therefore focus on the case where a < 0. Although the Time Dependent Ginzburg-Landau 
equation (/i — 0) is well known for phase separating systems [4] and as a generic amplitude 
equation [5] , we are not aware of previous work incorporating the convective term proportional 
to fi. This term breaks the variational character of the dynamics and the x — > —x symmetry. 
Since there exist two asymmetric kink solutions for any value of /i, we do not expect qualitative 
changes in the asymptotic dynamics of Ea. l|15|l when fi varies. Hence, we conjecture that the 
known result of logarithmic coarsening of Ea. (|15|) for /i = [13] extends to arbitrary /z. At 
short times, when /i is large enough, Ea. H15(l shares similarities with Burgers' equation, and 
its solution exhibits shocks. 

Let us now consider dynamics having a stable steady state h = 0, and a Lyapunov func- 
tional J- = J dx(j). The simplest dynamics which has this property (usually referred to as 
model A [4]), is: dth = —r{6J-/5h), where F is a function of h and its spatial derivatives. In 
presence of a small destabilizing non-equilibrium force /, we have: 

dth^-r^-^ + fK + o{f), (16) 

where K depends on h and its spatial and temporal derivatives. If J-' is translationally invari- 
ant, then is a function of dxh and its spatial derivatives only. Such a situation may be found 
in phase separation [4], where a gradient energy ~ J dx{dxh)'^ is used, or in the non-conserved 
dynamics of a thin film stabilized by surface tension. Linearizing Ea. H16|l . we then find that 
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e ^ f. Once again, an inspection of Ea. (|16|l proves that the dynamics is highly nonhnear (i.e. 
a = 0). To leading order in e, the evolution equation reads: (^) 



where P and Q are functions of h only. We shall notice here two unexpected and strongly 
restrictive properties of this equation: (i) it has the x ~x symmetry, although the dy- 
namics of the original problem does not necessarily have it. (ii) Although the front is not at 
equilibrium, Ea. H17(l exhibits a Lyapunov functional 



where i? is a function of h defined by the relation R' — Q/P. Indeed one can write dth = 
—PSU/Sh, and dtU < 0. Since Ea. (|17|l is variational, one can use the concepts of dynamical 
scaling developed for the study of phase transitions [4] to study the coarsening. One should 
once again rely on the structure of the steady-states in order to analyze the dynamics [19]. 

To conclude, we have presented two scenarios for the destabilization of a ID system: weakly 
nonlinear, and highly nonlinear dynamics. During weakly nonlinear dynamics, the evolution 
of the front morphology is described by an nonlinear expansion at small amplitudes. The 
amplitude then tends to zero as one gets closer to the threshold. During highly nonlinear 
dynamics, nonlincarities come into play only when the amplitude becomes finite. A small 
amplitude expansion is then not justified anymore. We have shown that the form of the 
evolution equation can still be found in this case. 
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